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We investigate the hydrodynamic interactions between microorganisms swimming at low Reynolds 
number. By considering simple model swimmers, and combining analytic and numerical approaches, 
we investigate the time-averaged flow field around a swimmer. At short distances the swimmer 
behaves like a pump. At large distances the velocity field depends on whether the swimming 
stroke is invariant under a combined time-reversal and parity transformation. We then consider 
two swimmers and find that the interaction between them consists of two parts; a dead term, 
independent of the motion of the second swimmer, which takes the expected dipolar form and a live 
term resulting from the simultaneous swimming action of both swimmers which does not. We argue 
that, in general, the latter dominates. The swimmer-swimmer interaction is a complicated function 
of their relative displacement, orientation and phase, leading to motion that can be attractive, 
repulsive or oscillatory. 

PACS numbers: 



Because of their size, the swimming motion of micro- 
scopic and mesoscopic organisms, such as bacteria and 
cells, corresponds to low Reynolds number. As pointed 
out by Purcell [l| , in this regime swimming mechanisms 
are very different from those operative at human length 
scales. In particular any swimming strategy must in- 
volve a cyclic and non time reversible motion. Examples 
include the helical motion of sperm flagellae ^2] , the syn- 
chronised beating of cilia in bacteria [3| , and the Dreyfus 
microswimmer [J| where a colloid is driven by a tail of 
magnetic particles actuated by an oscillating field. 

Interesting progress has been made in understanding 
the propulsion of individual swimmers and in 

developing optimal swimming strategies 8, 9, 10]. How- 
ever, much less is understood about the collective be- 
haviour of microswimmers. There is considerable evi- 
dence that groups of swimming cells undergo large den- 
sity fluctuations and form transient patterns (ill . fl^ . [isl . 
[l4 [TH . Indeed such collective behaviour appears to be a 
generic feature of so-called active systems, which produce 
their own energy [iGj . 

An important step towards understanding the be- 
haviour of many swimmers is to calculate the hydrody- 
namic interaction between two swimmers. Previous stud- 
ies have investigated this problem for a number of special 
cases. For example, Ishikawa and Hota [l3| performed 
experiments to find the interactions between two swim- 
ming Paramecium cells. Computational modelling of the 
hydrodynamic interaction between two flagella was re- 
ported in [l^[l^. Ishikawa et al. [20^ modelled swimmers 
as spheres with prescribed tangential velocities (known 
as squirmers). They calculated the far field interaction 
between two squirmers and matched the results onto a 
near field lubrication theory. An important disadvantage 
of this method is that it neglects any interaction resulting 
from the simultaneous motion of both swimmers. We find 
that these neglected terms are, in fact, very important for 
the dynamical behaviour of a wide class of swimmers. 

The aim of this article is to calculate two swimmer 



interactions, emphasising the importance of the details 
of the swimming stroke. We present both analytic and 
numerical results based on a minimal swimming model, 
first introduced by Najafi and Golestanian [21] . We first 
consider the average flow field around a swimmer. At 
short distances the swimmer behaves like a pump which 
moves forward by pushing fluid from in front to behind 
itself. At large distances the average flow field around a 
swimmer is naively expected to be dipolar (~ r~^). How- 
ever, we argue that for a wide class of swimmers, which 
are invariant under a combined time-reversal and par- 
ity (TP) transformation, the leading order term vanishes 
and the far flow field has a r~^ dependence. Examples of 
these swimmers include Purcell's three link swimmer fl | , 
the three sphere swimmer of Najafi and Golestanian [21 1, 
the snake swimmer and a helical filament. Indeed the 
term in the velocity field will remain important for many 
swimmers which are close to TP invariant. 

We then consider two swimmers, Alice and Bob, say. 
We find that the effect on Alice of Bob's swimming con- 
sists of two principal components; a dead term, which de- 
pends only on the flow field generated by Bob and which 
takes the expected dipolar form, and a live term resulting 
from the simultaneous swimming motion of both Alice 
and Bob which does not. Interestingly, it is the latter 
which dominates, meaning that interactions become not 
only a complicated function of relative displacement and 
orientation, but also of relative phase. We show that tra- 
jectories of the two swimmers can be attractive, repulsive 
or oscillatory. 

Locomotion is initiated by shape changes in the swim- 
mer, which induce forces acting on the surrounding fluid 
and thus create a flow. The fundamental solution of the 
Stokes' equations describing the response to a point force, 
f , is called a stokeslet and creates a fluid flow 

where y is the point of application of the force, r — |x— y 
and fl is the fluid viscosity. Since the Stokes' equations 
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FIG. 1; The cyclic motion of the hnear three sphere swim- 
mer [2^. The swimmer goes through the four steps (i) - (iv) 
and returns to its original shape, but displaced by a distance 
A to the right. 



are linear, the fluid motion associated with swimming 
may be given by a superposition of such stokeslet flows, 
which at large distances may be expanded in successively 
higher powers of and the far field flow approximated 
by the lowest order, non-zero term. Since the total force 
on the swimmer has to be zero the leading order term 
vanishes The next term, which varies as r~^, is called 
a force dipole and is g enerall y assumed to dominate the 
far field behaviour [H,!!!!!!. 

Indeed this is always true of the instantaneous fluid 
velocity. However, an intriguing feature of the average 
flow is that the force dipole term can be absent. To 
show this we employ a simple model swimmer which con- 
sists of three spheres of radius a joined by thin rods, or 
'arms'. It moves by shortening and extending these arms 
in a periodic and time irreversible manner, as depicted in 
Fig. [1] A unit vector n pointing in the direction of mo- 
tion, left to right in Fig. [U defines an orientation for the 
swimmer. The variable D gives the extended arm length 
and and denote the amplitudes of the strokes per- 
formed by the front and back arms respectively. Najafi 
and Golestanian [2l| obtained an analytic expression for 
the velocity of the swimmer for = using the Oseen 
tensor formulation of hydrodynamics, valid in the limit 
a/D < 1. 

Using a similar approach, we have calculated the aver- 
age flow field generated by the swimmer: 

+ [^^] ((1 - 3[n.f]^) ^ + (5[n.f]^ - 3) i^n^ll) } , 

(2) 

where T is the period of the swimming stroke. For 
swimmer is TP invariant and the dipolar 



term vanishes. (The amplitudes associated with each of 
the terms are given here only to lowest order in . 
However all terms in the coefficient of are antisym- 
metric in and hence vanish in the TP invariant 
limit.) 

The reason TP invariant swimmers do not have any 
dipole contribution can be seen as follows. If a general 
swimming motion is invariant under simultaneous time 
reversal and parity transformations then the average flow 
field must also have this property. The dipole term, how- 
ever, changes sign under a PT transformation and so the 
coefficient multiplying it must be zero. (This can be ver- 
ified for the first term in Eq. ([2]).) 

When the front and back arms have different ampli- 
tudes the TP invariance is broken and both and 
contributions appear in the averaged far field flow in 
Eq. ([2]). A comparison of the magnitude of the two con- 
tributions shows that the dipolar term only dominates 
for distances 
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which diverges when the amplitudes of the two arms be- 
come equal, as required by the TP invariance. 

The velocity field generated by this swimmer, averaged 
over a swimming cycle, is shown in Fig. [2l It was ob- 
tained using a numerical integration of the Oseen tensor 
equation, subject to the constraint that the swimmer is 
force and torque-free. Details of this method are given in 
[2H . Since the numerical technique does not exploit a far 
field expansion of the stokeslets it enables the flow near 
the swimmer to be determined as well, as illustrated in 
the inset. Close to the middle sphere the fluid is clearly 
seen to move to the left; essentially the swimmer acts 
like a pump which moves fluid from in front of it to be- 
hind it, leading to the swimmer moving to the right. The 
analytic far field expansion Eq. ^ accurately describes 
(to within 10%) the numerical flow for distances further 
than around three swimmer lengths. The dashed circle 
demarcates the estimated distance, Eq. ([3]), at which the 
far field flow crosses over from r^^ to dipolar scaling. 

It is interesting that two separate cases arise for the 
cross-over to a dipolar far field flow. This is because 
the term is antisymmetric under <-> while the 
term is symmetric, giving distinct flow fields for the 
two cases < (Fig.[2i;a)) and > (Fig.^b)). 
As we shall show, this distinction results in significantly 
different long time behaviour of interacting swimmers. 

We now consider how the trajectory of a second nearby 
swimmer (Alice) will be altered by this flow field. The 
interaction will involve Alice both being advected by, and 
rotating in, the flow field generated by Bob. We find that 
there are two contributions to the interaction. The first, 
which we refer to as the dead term, is independent of 
the motion of Alice. This contribution may be calculated 
using the average flow field for a single swimmer, Eq. 
which to leading order gives 



](^){3[nB.rp-l}, (4) 



3 




FIG. 2: The flow field around a swimmer averaged over one 
swimming cycle, (a) < with the inset showing the de- 
tail near the swimmer, and (b) > i^^. In both cases the 
direction of swimming is from left to right and the flow is 
axisymmetric about the swimmer. The magnitude of veloci- 
ties are plotted on a logarithmic scale to enable both the near 
sphere and far field behaviour to be observed simultaneously. 



The parameters used were a = O.ID and — 0.2D, 
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X |[nA.r] + 2[nA.nB][nB.f] - 5[nA.f][nB.f]2| , (5) 

where subscripts A, B label Alice and Bob, and r is the 
position vector of Alice relative to Bob. Axa repre- 
sents the translation vector, and AO a the rotation vector 
(whose magnitude and direction define the rotation angle 
and axis respectively), induced in Alice as a result of the 
hydrodynamic interaction with Bob. Although the rota- 
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FIG. 3: The function <1? describing the effect of the rela- 
tive phase of the swimmers on the interactions. The pa- 
rameters used were: black line = = 0.3-D, red line, 
= 0.473,5^ = 0.273, and blue line, = 0.2D, = 0.4D. 



tion is higher order in it should not be considered 
less important. This is because small changes in angle 
can give rise to large positional displacements in the long 
time limit. 

In addition, there is a live term which describes the 
more interesting interaction arising as a consequence of 
both organisms trying to swim. This term encodes all of 
the information about how the interactions depend on the 
relative phase of the swimmers. Using the Oseen tensor 
to describe the hydrodynamics we have determined the 
leading order contribution to the live interaction, finding 
a net swimming stroke translation and rotation given by 
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-5[nB.f]2-20[nA.nB][n^.f][nB.r]+35[nA.f]'[nB.f]'}. 

(7) 

The relative phase of the swimmers enters through the 
function $. The precise form of <i> depends on the amount 
of time taken for each of the four steps in Fig. [T] One 
choice is to assume that arm lengths vary at a constant 
rate w, so the time for steps (i) and (iii) is /w, and for 
steps (ii) and (iv) is ^^/w. In this case $ is plotted for 
a selection of parameter values in Fig. [31 Other choices, 
for instance requiring that each of the four steps takes 
equal time, produce very similar results. 

When the two contributions to the interactions are 
compared we can see that the live contribution domi- 
nates for all distances 
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This behaves similarly to the average fluid flow in that it 
diverges as the swimmer becomes TP invariant. However, 
there is a crucial difference for non TP invariant swim- 
mers. The additional factor of D/a is large, and therefore 
ensures that for all distances over which the interactions 
are strong enough to be significant, the dominant contri- 
bution is from the live term, despite the fact that it is 
higher order in . Thus, in general, even though the 
average flow field may be dipolar, the interactions are 
not. 

We turn now to an investigation of the long time be- 
haviour resulting from these interactions by considering 
how the trajectories of two identical swimmers depend 
on their relative position, orientation and phase. We fo- 
cus only on the most interesting case when the swimmers 
follow almost parallel trajectories since this allows time 
for the weak interactions to have an appreciable ef- 
fect. We consider first symmetric swimmers (^^ — 
swimming in phase, before looking at how the behaviour 
changes when these constraints are relaxed. The results 
are summarised in Fig.UI^a). Alice initially lies at the ori- 
gin (as shown in the figure) and Bob is aligned parallel 
and displaced to the corresponding point in the diagram. 
The figure is divided into different regions, labelled A, 
R, O, and P, corresponding to the four different types 
of long term behaviour that were observed; attractive, 
repulsive, oscillatory and parallel trajectories. 

Region A represents attractive behaviour, an example 
of which is illustrated in Fig. [H^i). Initially, Alice lies 
at the origin and Bob is displaced hD above her. The 
two curves show the centre of mass motion of the two 
swimmers in time. For accuracy, the trajectories were 
obtained numerically 25| , although the results which fol- 
low from iterating the analytic formulae, Eqs. ([4])-(l7]), 
are in close agreement. All scales are in units of D, 
the swimmer arm length. The parameters chosen were 
a = Q.ID and ^ = 0.31?, giving a displacement per cy- 
cle of A « (7/12)aCV-D^ = 0.005L». Thus, since the 
swimmers move a distance ^ 25D before colliding, this 
simulation lasted for 5000 swimming cycles, illustrat- 
ing that the swimming trajectories evolve on much longer 
times scales than that of the individual swimming stroke. 
The lines terminate when the swimmers get too close, as 
here the Oseen tensor approximation breaks down. 

Region R contains swimmers which repel each other, 
as shown in Fig. IH^ii) . The swimmers each rotate through 
a small angle of ~7°, with the rotation largely complete 
after ^--^3000 swimming cycles. 

Fig. m^iii) illustrates the oscillatory trajectories ob- 
served in region O. These oscillations take place over the 
course of 10000 swimming cycles, although this period 
varies strongly with the separation of the two swimmers. 
It is notable that the trajectories exhibit a slow drift su- 
perimposed onto the oscillations. The oscillations origi- 
nate in the rotational part of the live interaction, while 
the superimposed drift is a consequence of the transla- 
tional part. This highlights the significance of rotational 
interactions, despite their being higher order in r~^. 
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FIG. 4: Long time behaviour of two swimmers, Alice and 
Bob, (a) in phase, (b) tt out of phase, (c) with Alice phase 
lagging by 7r/2, and (d) in phase but with > (i)- 
(vi) illustrate a selection of trajectories representing attractive 
(A), repulsive (R), oscillatory (O), parallel (P), expanding 
oscillatory (OE), and contracting oscillatory (OC) behaviour, 
respectively. 
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Fig-HI^b) shows the diagram corresponding to Fig.[4l^a) 
but with the two swimmers tt out of phase. The trajec- 
tories change because the amphtude of the Hve interac- 
tion changes sign. The attractive region is replaced by a 
repulsive one. The oscillatory regime now produces os- 
cillating parallel paths, an example of which is shown in 
Fig. SJiv). Here the two swimmers follow parallel paths 
inclined at a small angle to the direction in which they 
would move in the absence of interactions. Both trajec- 
tories exhibit small amplitude, long period oscillations 
about the swimming direction; the oscillations of the two 
swimmers are exactly tt out of phase, and the amplitude 
of the oscillations vanishes when the angle between them 
is ~30°, corresponding to one of the zeroes of A0^^°. 

Fig. \M,c) shows the intermediate case when Alice has 
a 7r/2 phase lag with respect to Bob. This case is note- 
worthy as the amplitude for the lowest order live part of 
the interaction goes to zero. Nonetheless, it is clear that 
the live terms still dominate since the behaviour is dif- 
ferent for swimmers placed in front or behind the central 
swimmer. This positional change is equivalent to chang- 
ing the relative phase between the swimmers from 7r/2 
to 37r/2 and thus can only be important if the live terms 
are dominating the interaction. 

The oscillations of Fig.lD^iii) are perfectly periodic and 
repeat indefinitely. This is a direct consequence of the 
TP invariance of the swimmers when the amplitudes of 
the two arms are equal. When this invariance is broken, 
the oscillating states are no longer indefinitely stable. In 
particular. Fig. IDJv) shows the effect of making the front 
swimming amplitude, smaller and the back ampli- 
tude, , larger. Now the amplitude of the oscillations 
becomes successively larger until they can no longer be 
sustained and, finally, the two swimmers follow diverg- 
ing trajectories. The opposite case, when > is 
shown in Fig. IH^vi). The amplitude of the oscillations 
decays and the swimmers move closer together, with the 
lines terminating when the Oseen tensor approximation 
breaks down as a result of the swimmers getting too close. 
This behaviour is of particular interest as it suggests that 
suspensions of swimmers driven in phase might sponta- 
neously start to form chains. 

In this article, we have considered the collective dy- 



namics arising from hydrodynamic interactions between 
two model microswimmers. These interactions are not 
dominated by the leading non-zero term in a far field ex- 
pansion of stokeslets, but rather by a higher order live 
term describing the simultaneous motion of both swim- 
mers while they swim. Consequently, the long time be- 
haviour is sensitive to the relative phase of the two swim- 
mers and may be tuned by making small alterations to 
the details of the swimming stroke. 

Current experimental artificial microswimmers rely on 
external actuation in order to move 4]. Therefore, situa- 
tions in which two interacting swimmers are kept strictly 
in phase could be testable against the trajectories found 
in Fig. 2) In biological systems, the individual swim- 
mers move autonomously and the relative phases be- 
tween them would be expected to change in time. Given 
the sensitivity of hydrodynamic interactions to relative 
phase, such changes, perhaps leading to phase-locking, 
may have a significant impact on collective swimmer mo- 
tion and understanding this effect will be an important 
next step towards biological realism. 

The behaviour of two swimmers provides only a flavour 
of the collective dynamics of larger groups of microorgan- 
isms. Thus, although our analysis has shown that attrac- 
tive hydrodynamic interactions exist between two swim- 
mers, this fact alone is insufficient to claim clustering or 
swarming in larger groups of microorganisms. Similarly, 
the pairwise repulsion of two swimmers does not imply 
that larger groups will disperse. As such, there is much 
that still needs to be done to determine the collective 
behaviour of a large number of hydrodynamically inter- 
acting swimmers. This study has provided a framework 
in which this might be done and hinted at the richness 
which is likely to be found. 
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